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1 Introduction 



In this paper we construct the equilibrium Glauber and Kawasaki dynamics on the 
discrete particle systems such that certain determinantal point processes are invariant 
^ I under the Markov processes. 

5^ I The determinantal point processes, or fermion point processes, are point processes 

whose correlation functions are given by determinants of kernel operators. It was 
invented by Macchi [T3] and then has been extensively investigated by many people. 
It appears in many fields in mathematics and physics, for example, in random matrix 
theory and in fermion particle systems. For the details we refer to \10\ [20l [2T| [23] and 
references therein. 

The construction of equilibrium dynamics is one of interesting subjects for the de- 
terminantal point processes. One approach is to construct the diffusion processes via 
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Dirichlet forms |12j . During the last two decades there have been many works on the 
construction of equihbrium diffusion processes for Gibbs measures of interacting parti- 
cle systems. We notice that the interacting Brownian particles of logarithmic potential, 
which is related to the determinantal point process of sine kernel and is called Dyson's 
model, was studied by Spohn [25j. Some related works to Dyson's model were also 
done recently in [SKTB]. The diffusion process via Dirichlet form for the determinantal 
point processes in the non-percolating regime, or equivalently in the high tempera- 
ture or small activity regime, was constructed by the second named author [25] . The 
other approach, in particular for particle systems, is to construct the particle birth 
and death processes, so called Glauber dynamics, and the particle jump processes, 
called Kawasaki dynamics. The general interacting particle systems in discrete model 
were developed in detail by Liggett [9]. The Glauber and Kawasaki dynamics for the 
continuum models have been investigated in the literature [31 HI El [3 El E] ■ Among 
these, Kondratiev et al. established the Dirichlet form approach for the construction 
of the equilibrium Glauber and Kawasaki dynamics for continuum systems so that the 
Gibbs measures for the system are invariant under the Markov processes [3 [8] . For 
instance, the standard superstable interaction of Ruelle [18] falls into their regime of 
applications. Lytvynov and Ohlerich applied the methods developed in [71 [8] to con- 
struct the equilibrium Glauber and Kawasaki dynamics that leave invariant certain 
determinantal point processes in continuum model [TT]. The Glauber dynamics for 
discrete determinantal point processes was studied by Shirai and the second named 
author [22j . This paper can be regarded as a continuation of the work in [22]. We 
emphasize that the Markov processes constructed in this paper are Feller processes. 
That is, the semigroup acting on a continuous function gives rise to another contin- 
uous function. But the Markov processes constructed via Dirichlet forms are Hunt 
processes and the semigroups act on L^-functions. 

Let us briefly sketch the contents of this paper. We consider the infinite particle 
system, particles living on a discrete set, say E. We consider the determinantal point 
process /i on the configuration space X. The point process /x has a defining kernel 
operator K of the type K := A{I + A)~^, where A is a positive definite, bounded 
linear operator on P{E) that satisfies some hypothesis (H) (see Section 2). Under 
the hypothesis (H), /x is known to be Gibbsian [26]. The generators of Glauber and 
Kawasaki dynamics have the form (see subsection 2.1 for the details): 

and 
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where d{x;^) and b{y; ^) are death and bnth rates for Glauber dynamics, and c{x, y; ^) 
are the jump rates for Kawasaki dynamics. In Subsection 13. H following Liggett [9], 
we introduce the basic strategy to show the existence and ergodicity of the dynamics. 
Next for point processes with Papangelou intensities, we will give necessary conditions 
for the rate functions to satisfy the detailed balance condition (Theorem 13. 3p . Then 
we apply these ideas to our model. Given a DPP /i for a kernel operator which satisfies 
certain hypothesis, we give some concrete formulas for the rate functions (Proposition 
13. 5p and present the conditions for existence in the language of Papangelou intensities 
(Proposition 13. 7p . To guarantee that the conditions for the existence and ergodicity 
are satisfied, we need further stronger conditions introduced in Assumption (A) in 
Section 4. Under the assumption (A), we finally construct the dynamics of our purpose 
(Theorem l4.2p . In the Appendix, we provide with a proof of a technical lemma, Lemma 
14. H which is worth to be noticed in itself. 

Finally, comparing to [22], we would like to mention that in this paper the Kawasaki 
dynamics is included and the concept of Papangelou intensities is used in a crucial 
way to construct the dynamics. After denoting the detailed balance condition by Pa- 
pangelou intensities, it is possible to choose the flip or jump rates for the equilibrium 
dynamics in many ways. In [22], we dealt with just one choice among them (see Re- 
mark [3i6|). However, at the moment, we have to assume seemingly almost the same 
conditions as in [22]. Our future studies are addressed to 

- the construction of the Glauber and Kawasaki dynamics under weaker condi- 
tions, at least under the hypothesis (H); 

- the investigation of the spectral gap, or log-Sobolev inequalities for the genera- 
tors of the Markov processes. 

2 Preliminaries 

In this section we briefly recall the definition of Glauber and Kawasaki dynamics for 
spin systems, or equivalently lattice gases. Then we introduce the determinantal point 
process in discrete spaces and their Gibbsianness. 

2.1 Glauber and Kawasaki Dynamics for Lattice Gases 

Let E be any countable set. We have in mind the system of the lattice space with 
E = Z"^, the (i-dimensional integer lattice, but we deal with arbitrary discrete set E. 
Let X be the set of all subsets C E, called the configurations. For any subset A C E 
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and ^ € we denote by the restriction of ^ to A: 

U-=CnA. (2.1) 

From now on if a subset A C is finite, we denote it by A CC i?. For each A C -E, let 
Ta be the a-algebra on X generated by the functions ^ i-> |^a|, A CC A, where |^a| 
means the cardinahty of the set |^a|- Thus we get an increasing system of cr-algebras 
{-^a}acc-B and we let T := Te- Any probability measure ^ on {X^F) is called a 
point process. We notice that the cr-algebra T on X can be understood as a Borel 
cj-algebra by a natural identification of X with {0, 1}^ equipped with the product 
topology. This identification is taken for granted in this paper, and consequently we 
will also consider ^ as a topological space. Notice in particular that X \s a compact 
space. 

Let C{X) denote the set of all continuous functions on X equipped with the sup- 
norm II • ||. For / G C{X) and x ^ E, let 

A/(x) := sup{|/(xe) - fm : e e ^}, (2.2) 

where we have used a shorthand notation := {x}U^. We define a subset of "smooth" 
functions D{X) [9! by 

DiX) := {/ € CiX) : |||/||| := ^ A;(x) < oo}. (2.3) 

By the Stone- Weierstrass theorem, we easily see that D{X) is dense in C{X). 

The generators for Glauber and Kawasaki dynamics are defined as follows. We 
first consider Glauber dynamics. For each x (z E, let b{x; ^) and d{x; ^) be nonnega- 
tive continuous functions on X. They are called birth and death rates, respectively. 
Namely, given a configuration ^ E ^, a particle is born at site x ^ ^ in a rate 6(x;^), 
and among the particles ^ a particle at x € ^ dies out in a rate d{x; S,). For each x G 
we define 

Cx := supmax{6(x; ^), d(x; x^)}, (2-4) 

i^x 

and define 

c(^):=supc^ (2.5) 

Next for Kawasaki dynamics, we need the particle jump rates. For each x ^ y (z E, 
let c(x, y; ^) be a continuous function on X whose values are defined to be zero unless 
X G and y ^ ^. It is the rate for a particle at x G ^ to jump to the empty site y G E\S^. 
We will also need to control the particle jump rates. For each pair x ^ y £ E, define 

C{x,y} ■= sup max{c(x,y;xO,c(y,x;y^)}; (2.6) 
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:=supVc|,,j,}. (2.7) 

We will assume that c*^"^ < oo for [J = G or K. Under this assumption, the Markov 
pregenerators for Glauber and Kawasaki dynamics are defined as follows. For / G 
D{X), we define for Glauber dynamics 

and for Kawasaki dynamics 

i^^V(e):= E (2-9) 

Here we have also used a short notation (,\x for ^ \ {x}. For a definition and proof 
that L^'^\ [j = G or K, becomes a Markov pregenerator, we refer to [9l Definition 2.1 
and Proposition 3.2]. 

In this paper we will investigate the conditions on the rates 6(x;^), d{x;^), and 
c{x,y;(^) so that not only (the closure of) L^^\ {j = G or K, becomes a Markov gen- 
erator on C{X) but also it leaves invariant certain point process /.t on {X,J^). We 
are concerned exclusively with determinantal point proccesses on E, which are briefly 
introduced in the next subsection. 



2.2 Determinantal Point Processes 

Determinantal point processes (DPP's hereafter), or fermion point processes, are the 
probability measures on the configuration space of particles. The particles can stay 
either on continuum spaces or on discrete sets. The correlation functions of DPP's 
are given by determinants of a priori given kernel operator as shown in Theorem 12.11 
below. Typically, they have a fermionic nature, namely, the energy increases when 
new particles add into a given configuration of particles. For the basic theories of 
DPP's we refer to references \10 \ \13 \ [2T \ I23j. Here we follow the reference |21] for the 
introduction. 

Let E be the countable set in the previous subsection and let K he a Hermitian 
positive definite, bounded linear operator on the Hilbert space -Ho := P{E), the set of 
square summable functions (sequences) on E equipped with the usual inner product: 

{f,9)o:=Y.^9{x), f,g£no. (2.10) 

The induced norm is denoted by || • ||o. The following is an existence theorem for 
DPP's, which we present in the form stated in |21j . 
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Theorem 2.1 Let E be a countable set and K a Hermitian bounded operator on Hq. 
Assume that < K < I . Then, there exists a unique probability measure fi on 
such that for any finite set {xi, . . . , Xn} C E, the n-th correlation function is given as 
follows: 



p'i^'{xi,...,Xn) := G : C :5 {xi, . . . ,x„}}) = det{K{xi,Xj))i<ij<n. (2.11) 



Next we discuss the density functions for DPP's. For each subset A C E, let Pa 
denote the projection operator on T-Lq onto the subspace /^(A) and let K\ := P\KP\ 
denote the restriction of K onto the projection space. For each A CC -E, assuming at 
the moment that I\ — K\ is invertible, we define 



The local marginals of the DPP fi corresponding to the operator K are given by the 
formula: for each A CC E and fixed ^ € A", 



/^a(Ca) := K{C ■ Ca = Ca}) = det(/A - Ka) det{A[f^]{x,y))^^y(z^j^, (2.13) 



where A[A](x,y), x,y € A, denotes the matrix components of ^[aj- We remark that 
the r.h.s. of (j2.13p can be given a meaning even when I\ — K\ is not invertible j21U23j . 

2.3 Reproducing Kernel Hilbert Spaces, Papangelou Intensities and 
Gibbsianness of DPP's 

In this subsection we briefly discuss the Gibbsianness of DPP's. To show the Gibb- 
sianness of a point process fj, is equivalent to show the existence of global Papangelou 
intensities of fi \14:\ [T5\ \T9\ [2T] . Since the Papangelou intensities are the key ingredients 
for the construction of equilibrium dynamics we review it for our model from |26j . We 
start by introducing a dual pair of reproducing kernel Hilbert spaces pQ. 

Let A be a positive definite, bounded linear operator on T-Iq = l'^{E) (the operator 
will define a DPP via the operator K := A{I + A)~^). We assume that Ker^ = {0}, 
so Ran^ is dense in Hq. We introduce two new inner products (•,•)- (')')+) 
respectively on T-Lq and Ran A as follows. First on T-Lq, define 



^[A] := Ka{Ia - Ka) 



1 



(2.12) 



if, 9) 



{f,Ag)o, f,geno, 



(2.14) 



and on RanA define 



(/,5)+ := {f,A-' 



g)o, f,g£RanA. 



(2.15) 
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The induced norms will be denoted by || • ||„and || • ||_|_, respectively. Let and T-L^ 
be the completions of T-Lq and Ran^ by the respective norms || • ||_and || • || + . We then 
get the following rigging of Hilbert spaces: 



We let B := {ex}x£E be the usual basis of Tio, i.e., Cx is a function on E whose 
value at X is 1 and the values at other sites are all zero. Let A{x,y),x,y G E, be the 
representation of A w.r.t. B. Then we notice that is a reproducing kernel Hilbert 
space (RKHS hereafter) with reproducing kernel (RK shortly) A(x,y). That is, ■H+ 
is a linear space of functions on E with the following defining conditions: 

(i) For every x ^ E, the function A{-,x) belongs to %+] 

(ii) The reproducing property: for every x E and g € 7^+, g{x) = {A{-,x), g)^. 

We want H- to be also a RKHS (it is not the case in general), so we assume the 
hypothesis below throughout this paper: 

Hypothesis (H): We suppose that is functionally completed [1], i.e., any vector 
of Ti^ can be represented as a function on E. 

For sufficient conditions for the hypothesis (H) we refer to \26\ Appendix]. Under the 
hypothesis (H), becomes also a RKHS with RK, say B{x,y),x,y G E [Tl [26]. 
Informally speaking, B{x,y) is the matrix representation of A"^, though it is not 
of bounded operator in general. In particular, under (H) we have G 71+ for all 
X £ E The fiip or jump rates of our dynamics will be defined via the Papangelou 
intensities (defined below) of the DPP /x corresponding to the operator K := A{I + 
A)~^, which are in turn the squared norms of the vector-projections in the Hilbert 
space H- [26]. In particular, the interdependencies of rates, which must be controlled 
for the construction of the dynamics, are represented by the inner products in the 
restricted Hilbert spaces of H- (see Proposition 13. 8p . 

As mentioned above, we are concerned with DPP's corresponding to the operators 
K := A{I + A)~^ . In order to get the Papangelou intensities, we need a variational 
principle [26]. For each A CC E, let 




Let xq £ E he a fixed point and let E = {xq} U RiU R2 be any partition of E (one 
of Ri and R2 might be the empty set). For each A CC E, define 



(2.16) 




Fioc,A '■= the linear space spanned by {ex ■ x € A}. 



(2.17) 



OA : = 




/II- and /3a: 




(2.18) 
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Obviously, {aA}Acc-B arid {/3a}acc-E are decreasing nets of nonnegative numbers. 
Consequently we define 

a:=limQA and B := lim B\. (2.19) 

The variational principle says that, under the hypothesis (H), no matter how we take 
a partition E = {xo}Ui?iUi?2, the product of a and /? is equal to 1 (see [26j Theorem 
2.4] and also [M]): 

a(3 = 1. (2.20) 

The relation (|2.20p guaranties, on the other hand, the existence of global Papangelou 
intensities. Let fi be the DPP corresponding to the operator K := A{I + A)~^ , where 
A satisfies the hypothesis (H). Recall that the local Papangelou intensities of fi is 
defined as follows: for each x G E, x G A CC E, and x ^ ^ G X, 

The global Papangelou intensities are the limits 

a[x;i) := lima[A](x;^A), (2.22) 
whenever the limit exists. The following theorem was proved in [26t Theorem 2.6]: 

Theorem 2.2 Let the operator A satisfy the hypothesis (H) and let be the DPP 

corresponding to the operator K = A{I + A)""^ . Then for all x G E and x ^ G X , 
the Papangelou intensity a{x; ^) exists and it is equal to the the number a in ()2.19p 
obtained by replacing xq and Ri in ()2.18p by x and respectively. 

Remark 2.3 By (fZTSD and (f2:2n) - (f2:22D we see that 

det Ai\](xfA,xfA) 
a x;0 = lim . ' , , ^.23 

AtE det A[a](Ca,^a) 

here ^[a](Ca;Ca) is the matrix (^[A](2^)2/))x,j/eCA- What we have shown in the above 
theorem is that it is equal to the limit 

. detA{x^A,x^A) 
a(x:t) = hm — — — — . (2.24) 

We will denote the dual relation (j2.20p as 

a(x;0/3(x;0 = 1- (2-25) 
We notice also that for all x E and x ^ ^ X , 

a{x;^) < A{x,x) <\\A\\. (2.26) 
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3 Construction of Glauber and Kawasaki Dynamics 

In this section we will construct the Glauber and Kawasaki dynamics for DPP's. We 
begin by briefly introducing the general strategy for the existence and ergodicity of 
the dynamics following [9]. 

3.1 Existence and Ergodicity 

Throughout this subsection we assume that the flip rates d{x;(,), and jump 

rates c{x, y; ^) are continuous functions for ^ € <-f and satisfy the boundedness condi- 
tions c^*^ < oo, 7^ = G or K. Under these conditions the operators L^*), # = G or K, 
in (j2.8p - (j2.9p are Markov pregenerators [H Proposition 3.2, Chapter I]. 

The Markov pregenerators are closable [HI Proposition 2.5, Chapter I], but in order 
that the closures to be Markov generators, we need to control the interdependencies of 
the rate functions. For Glauber dynamics, following [9], we deflne for each x ^ u E, 

-f^^\x,u):= sup{\b{x;^) -b{x;u^)\ + \d{x;xu£,) -d{x;xC)\}. (3.1) 
And we deflne 

M(G) := sup V7(^)(x,'u). (3.2) 

We will need also the lower bound of the flip rates deflned by 

e(^) := inf inf {6(2;; C) + d{x] x^)}. (3.3) 

The interdependencies for Kawasaki dynamics will be given in the following way. First 
deflne for each x ^ u ^ E 

-f^^\x,u) ■.= ^ sup {\c{x,y-xS) - c{x,y]xuC)\, \c{y,x;y^) - c{y,x;yu^)\}. (3.4) 
Then we deflne 

m(k) :=supy7(K) {x,u). (3.5) 

u^x 

The lower bound of the rates for Kawasaki dynamics is given by 

e(^) := inf inf { V c(x, y; + E y^)}" (^.6) 

x^i x^yi 

Deflning 7^^) (x, rc) = 0, we let r(#) = (:i^*\x, u))^,ueE for # = G or K. The following 
theorem is proven by Liggett (9l Theorem 3.9, Chapter I]: 

Theorem 3.1 For # = G or K, assume that c^*) < oo and M(#) < 00. Then we 
have 
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(a) 



The closure U*) of L^*^ is a Markov generator of a Markov semigroup 



-{#) 



(b) 



D{X) is a core for 



(c) 



For f (^D{X) 



exp(tr(#))A/. 



(d) 



If f e D{X), then T^^^'f G D{X) for all t > and 



r(#Vlll<exp[(M(#)-£(#))t] Ill/Ill. 



Let us now introduce the concept of ergodicity. Recall that a point process on 
is said to be invariant for the Markov process with semigroup {Tt,t > 0} if 



for all / G C{X) and i > 0. A Markov process with semigroup {Tt,t > 0} is said to 
be ergodic [S] if 

(i) there is a unique invariant measure, say z^o; 

(ii) liuit^oo^'Tt = for all probability measure v on {X,J^). 

Theorem 3.2 f'JPl Theorem 4-1, Chapter I]) Suppose that the same conditions as in 
Theorem lg.il hold. In addition, if M^*) < e^*) then the process is ergodic. Further- 
more, for g G D{X), 



where v'^> is the unique invariant measure. 
3.2 Detailed Balance Condition 

In this subsection we discuss the detailed balance conditions for Glauber and Kawasaki 
dynamics. Let /U be a point process and suppose that ^ has Papangelou intensities 
q(x;^). By the general theory, this implies that /x satisfies the DLR-conditions, which 
read as follows: for all bounded measurable functions F : A" — )• M and A CC i?. 



where ^a(C) = E^acA "(Ca; CaO and a(CA; CaO = nl=i "(a^i; a^i • • • Xi^i^t^c). We refer 
to [m [15] for more details. 






(3.7) 
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The detailed balance condition for Glauber and Kawasaki dynamics w.r.t. fj, means 
that the pregenerators L^'^\ [| = G or K, are symmetric: 

J = I fiidOL'^^^fiOoiO, f\9 e D{X). (3.8) 

From the DLR-conditions (|3.7p it is not hard to get equivalent conditions for the 
detailed balance. For simplicity, we assume that the Papangelou intensities a{x; 
are positive. 

Theorem 3.3 Let L^^\ ^ = G or K, be the Markovian pregenerators for Glauber 
and Kawasaki dynamics given in (|2.8p and (|2.9p . In order that the detailed balance 
condition (j3.8p is satisfied, it is necessary and sufficient that the rate functions satisfy 

(a) for Glauber dynamics: for allx^E,x^S,£X, 

b{x;(,) = a{x;^)d{x;x(,); (3-9) 

(b) for Kawasaki dynamics: \/x ^ y £ E , and € <-f with x,y ^ ^, 

a{x] C)c{x, y; xO = a{y; ^)c{y, x; y(). (3.10) 

3.3 Glauber and Kawasaki Dynamics for DPP's 

In this subsection, we will concretely define the flip and jump rates for Glauber and 
Kawasaki dynamics for DPP's, and then discuss the existence and ergodicity condi- 
tions for the dynamics. 

From now on we fix an operator A on T-Lq that satisfies the hypothesis (H). Thereby 
we also fix a DPP /x with kernel operator K := A{I + A)~^. We will construct 
Fellerian Markov generators for Glauber and Kawasaki dynamics on C{X) that leave 
invariant the above DPP /x. Recall the notations in subsection 2.1. In order that the 
pregenerators (which are not yet explicitly defined) L^*-^) and L^^^ in (j2.8p and (j2.9p 
would define Markovian generators for Feller processes, it is needed that the functions 
b{x;(,), (i(x;^), and c(x,y;^) are continuous functions for ^ £ X. These flip, or jump 
rates will be given via the Papangelou intensity function a(a;;^). So, we need the 
following 

Lemma 3.4 For each x £ E, the Papangelou intensity a{x;(,) of fi is a continuous 
function of ^ € X. 
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Proof. Recall the definitions a a and /?a in (j2.18p . For each x ^ E and ^ ^ X with 
X ^ ^, we consider the partition E = {xjU^U^, where ^ = For each K dd E 

we define local functions ai^{x;S,\) and /3a(x;^a) by the formula (j2.18p replacing xq 
and i?i by x and ^, respectively. As local functions, clearly the functions ^ i— )• aA(x; ^a) 
and ^ I-)- f3\{x;£^\) are continuous on the set G Af : .^^x}. Now as decreasing limits 
of continuous functions, the Papangelou intensity a(x;^) and I3{x]^) are both upper 
semi-continuous functions on € A:* : C^x}. Since a(x;^) and /3(x;^) are reciprocal 
to each other, they are also lower semi-continuous. This proves the lemma. □ 

In the rest of the paper we will use the rate functions given in the following 
proposition, which we can easily prove. Recall the variational relation a(x; ^)/3(x; ^) = 
1 for all X G and X ^ C G A' in ([2:25]) . 

Proposition 3.5 Assume the hypothesis (H) and let be the DPP corresponding to 
the operator K := A{I + A)""^ . Then the following choices for flip rates b{x;^) and 
d{x; ^) for Glauber dynamics, and jump rates c(x, y; ^) for Kawasaki dynamics are 
uniformly bounded for G X and satisfy the detailed balance conditions w.r.t. fi: (we 
let x,y a) 

(a) for Glauber dynamics: 

:= ^MI- and d{x;xO := (3.11) 
^ l + a(x;0 l + /3(x;e) ^ ' 



(b) for Kawasaki dynamics: 

c(x, y; xi) = d{x, y)a{y- i)g{a{x\ C),a{y\ 0), (3.12) 

where d{x, y) is a symmetric weight function, and g : — )■ M+ is any symmet- 
ric, continuous function that makes c(x, y; x^) bounded. 

The most simplest example for the function g in the above is g{u,v) = 1. However, 
reflecting on the nature of the dynamics, including this case, we may choose for any 
< t < 1, 

5(n,.)=^^,(u,.):=(^^^-J^^) . (3.13) 
Then c(x,y;x^) becomes 
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Remark 3.6 In the equations ()3.1ip and ()3.14p . the terms (1 + /3(x; ^))~^ and ((1 + 
I3{x; C))(l + /3(y; to make the flip rates bounded. When a{x; ^) is uniformly 

(for ^) away from 0, we may drop these terms (taking b(x; ^) = a{x;^) and d{x; ^ = 1 
for Glauber dynamics). In [22], we have taken, in our terminology, b{x; ^) = l + a{x; ^) 
and d{x;x£,) = 1 + /3{x;(,) (see [22l eq. (1.9)]). 

Now the flip rates for the equilibrium Glauber and Kawasaki dynamics have been 
given via the Papangelou intensities a{x;(,), we would like to represent the condition 
M^*^ < oo in terms of a{x; ^). Let us assume that the weight function d{x, y) in (j3.12p 
satisfies 




(3.15) 



Proposition 3.7 Suppose that the flip rates b{x;(^) and d(x;^) and the jump rates 
c{x,y;x(^) are given by the equations p.lip - ()3.12p with the function g in (|3.12p given 
by p.lSp . and d{x,y) satisfying (j3.15p . Define 

Mj*^^ := sup sup (a(x; ^) — a(x; n^)) (3.16) 

and 

m['^ := SM^^^d{x,y) sup [{a{x;C) - a{x]uC}) + {a{y;C) - a{y]uC})] ■ 

(3.17) 

Then M(#) < aoM^"^^ with a uniform constant oq, and hence if M-[*^ < oo then 
< oo for ^ = G or K, and all the results in Theorem VJ.1\ hold. 

Proof. The proof immediately follows from the definitions. In particular, for the 
Kawasaki dynamics, we use the mean value theorem and the boundedness of the 
partial derivatives for the function gt{u,v) = {{1 + u)[l + v))~^ . □ 

(-U-\ 

In order to get the boundedness of M{ , # = G or K, we have to know the 
quantities a(x; .^) — a(x; uE,) in (j3.16p - (|3.17p more in detail. For that purpose we briefly 
introduce the restriction theory of reproducing kernel Hilbert spaces [1]. Let T-L be any 
RKHS (on E) with RK C{x,y). Let R C E he any (finite or infinite) subset of E, 
and let Cfi{x,y), x,y £ R, denote the restriction of C to the set R. It was shown by 
Aronszajn that Cr{x, y) is the RK of the RKHS, call it T-Lr^Cri on the set R consisting 
of all restrictions of H to the set R [T]. We let || • ||c and || • \\r^Cr denote the norms 
on T-L and T-Lr^Cr, respectively. Then || • \\r,Cr is given by for each / G 'Hr^Cr-, 



WfU^CR ■■= mf{|l/|b : fix) = fix) for all x € R}. (3.18) 



14 



Chae and Yoo 



Moreover, for each / € Hr^Cr^ there is a unique vector f Ti whose restriction to R 
is / and 

\\f\\R,Cn = ll/'llc (3.19) 

By the restriction theory for the RKHS's, the key terms a{x; ^) — a(x; n^) in the above 
proposition have the folfowing representations (see \'27\ Proposition 3.2]): 

Proposition 3.8 For any x ^ u ^ E and x,u ^ C E, we have 

a{x;^) -a(x;0 = | (e^;, e„)gc^s^, ^ ||e„||^,^^^^ . (3.20) 

In particular, in a formal way, we also have the representation: 

a{x;S,) -a{x;u(,) = \A{x,u) - ^(x,0^(C,C)-M(^, n)^ a{u;(,)-^. (3.21) 

In the next section we will discuss some sufficient conditions so that we can control 
(^{x] S,) ~ Oii^] u^) uniformly for ^. 

4 Examples 

In this section we discuss some examples for which the resulting Glauber and Kawasaki 
dynamics are Fellerian Markov processes on C{X) and leave invariant certain DPP's. 

(44:) 

In Proposition l3.7l we have seen that if we could control the interdependencies; < 
oo, then we are done. Unfortunately we couldn't do it under our hypothesis (H), so 
we impose further stronger conditions on the operator A. For any complex number 
z € C, we let \z\i := \Rez\ + |Imz|. 

Assumption (A): Let yl be a positive definite, bounded linear operator on T-Lq = 
l'^{E). We assume that there is a A > such that 

inf (a(x,x)- V|^(x,?/)|i) >A. 

x&E \ ^ — ^ / 

Any operator A that satisfies the Assumption (A) is said to be diagonally domi- 
nant. When one considers the convolution operators on it is not hard to see 
that there are many operators A that satisfy the assumption (A). For example, let 
C{x,y) = C{x — y) be the convolution operator on /^(Z'^) coming from the Fourier 
coefficients C(-) of a sufficiently smooth positive function p{t) on the torus T*^. Then 
the components C(x, y) decrease fast enough so that |C(x, y)\ < oo. We may then 
take A:= al + C for some positive numbers a. 
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Throughout this section we suppose that our operators A satisfy the assumption 
(A). Let us fix a constant g > so that 



q> q{A) := sup J2\Mx,y)\i- (4.1) 
We define a Q-matrix Q on E by 



. -\A(x, v)\i, X ^ y 

Q(x,y) = r (4-2) 

-Ey^xg\Mx,y)\l, X = y. 

Let n be a stochastic matrix on E defined by 

U:=Q + I. (4.3) 

The following lemma is an analogue of [22^ Lemma 4.1]. The proof follows by modifying 
the idea of the proof in [22] and we leave it at the appendix for readers' convenience. 

Lemma 4.1 Assume that the operator A satisfies the assumption (A). Then for any 
C E, the submatrix A{S^,^) is invertible and for any x,y € S,, 



m,0'\x,y)\<Mix,y) 



J, x = y, 



(4.4) 



where T := ^^^^ j^U 



With the help of the above lemma, we are able to state our main result. 

Theorem 4.2 Suppose that the operator A satisfies the assumption (A). Let us take 
the flip rates 6(a;;^), d(x;^), and c(x,y;^) as stated in Proposition \3. 7\ Then the 
hypotheses of Theorem \3.1\ are satisfied and hence all the statements of Theorem \3.1\ 
hold and the DPP fi corresponding to the kernel operator K := A{I+A)^^ is invariant 
under the Glauber and Kawasaki dynamics. Further, in addition to it, if q{A) in (j4.ip 
is sufficiently small, then the Markov processes are ergodic and the statements in 
TheoremWMhold. 



Proof. For the existence of Markov processes for Glauber and Kawasaki dynamics, by 
Proposition EZl it is enough to check that m[*\ # = G or K, defined in (l3T6]l - (l3Tn ). 
are finite. Under the assumption (A), ()3.2ip has a rigorous meaning. By Lemma [4.H A 
and any submatrix of A are boundedly invertible and we have ||^(C; C)~^ II ^ (-^+9) 
uniformly for ^ ^ E. Then we easily see that a{x;£) = {^A[x^,x^)~^{x,x)^ ^, and 
Lemma |4 . 1 1 gives the bound: 

0(2;; > A, uniformly for x € £^ and x ^ ^ C E. (4.5) 
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By using again the result of Lemma 14.11 in (j3.2ip we see that 

snj>\A{x,OA{^,0-'A{C,u)\ < sup ^ Kx,y)|M(y,z)|A(z,. 



< Yl \A{x,y)\M{y,z)\A{z,u)\. (4.6) 

y,z£E\{x,u} 



Hence we have by ^^2T\i and (fi3]) - (fO]) 

m}*^^ = sup sup {a{x;^) — a{x; u^)) 

< sup^(|A(x,n)|+ Yl \A{x,y)\M{y,z)\A{z,u)\)^ 



< ('-''^)M('-'^)<- 

(K) 

In a very similar way it follows that M\ < oo. 

The flip or jump rates are chosen so that the Markov processes satisfy the detailed 
balance conditions w.r.t. the DPP So, [i is reversible for both Glauber and Kawasaki 
dynamics, and hence [i is invariant under the dynamics (see O Propositions 5.2 and 
5.3, Chapter II]). 

In order to check the ergodicity we have to know the quantities e'^*^ # = G or K. 
By using the definition of the flip and jump rates and also by using the property (j4.5p 
it is not hard to check that 

e(G) = i and > A ^^^-i^, (4.8) 

where d\ is the constant in (j3.15p and is the operator norm of A. By Proposition 
13.71 (|4.7p . and (j4.8p . we see that if q is small enough, then M^*) < e(*\ and this 
completes the proof. □ 



A Appendix 

In this appendix we provide with the proof of Lemma 14.11 The main ingredients of 
the method are the Markov chain on the discrete set i?, Feynman-Kac formula, and 
the comparison of stochastic matrices of the Markov chains. The central idea was 
introduced in the proof of [221 Lemma 4.1]. 

Proof of Lemma \4-l\ We flrst assume that A is a real matrix. For each ^ C E, we 
define a Q-matrix on the set .^U^U {9}, where d is an abstract extra point playing 
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as a cemetery (when ^ = E, we just ignore it), and ^ is also an abstract set consisting 
of the elements of copies of we denote them as ^ := {x : x € ^}. We define 



Q^{x,y) = Q^ix,y) 

Q^{x,x) = Q^{x,x) 
Q^{x,x) = Q^{x,x) 

Q^{x,d) = Q^{x,d) 
Q^{d,x) = Q^{d,x) 



-A{x,y)-, x^y(^i, 

q 
q 

0, X G e, 

- V -\A{x,y)\, xee, 
a 

-\A{x,y)\, xet 
Qi:{d,d) :=0, xet 



(A.l) 



here a± := maxjita, 0} for each real number a. Given a function / € we extend 

it to a function / G /°°(^ U ^ U {d}) by 



fix) := <^ 



/(x), x = x 
— /(x), x = x 
0, X = d. 



(A.2) 



We consider the anti-symmetric subspace l'^{£,U {d}) of /°°(^ U ^ U {d}) defined 

by 

C(eueu{5}) := {/ G r{^U^U{d}) : /(x) = -/(x), /((?)= 0}. (A.3) 

We notice that can be regarded as an operator on U ^ U {d}). We define a 
function V^(x) on ^ U ^ U {5} by for x G ^, 

Vi:{x) = V^{x):=-(^-A{x,x)+ Yl \Mx,y)\), s^nd V^{d) = 0. (A.4) 



Now we consider the equation: 



(A.5) 



for h G By the definition of the Q-matrix Q^, the equation (lA.Sp is equivalent 

to 

(A.6) 



q 

Let us define a stochastic matrix on ^ U ^ U {(9} by 

:= + -^5u5u{a}- 



(A.7) 
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Let {xj:^^ : t > 0} be the Markov chain on U U {d} generated by Q^. Then by the 
Feynman-Kac formula, we have 

1 



fix) = -^iQ^ + nr'h{x) 



1 



^"/i(X«))exp V^{XP)dsjdt 



(Ai 



Let us now take h := 6y, the Dirac function at the point y. Then the solution f{x) in 
(jA.Sp is f{x) = A{^,^)~^{x,y), thus we have (notice that 6y = Sy — 5y) 



i^(e,e)-'(x,y)i<^E, 



5,MxP))exp( / V^{7,{X]^^))ds]dt 



where vr is a projection operator on ^ U ^ U {d} defined by 

Ti{x) = 7r(ir) := x and vr(9) := d. 



(A.9) 



(A.IO) 



We would like to estimate the r.h.s. of ()A.9p . For this, we introduce another Q-matrix 
on e U {d} by 



Q^{x,d) 



X ^y, x,y 



-^EzeE:z^x\Mx,z)\, x = y£C, 



(A.ll) 



We also define a stochastic matrix 11^ on ^ U {d} by 

:= + I^^^oy. (A.12) 

Then we notice that the probability law of the chain {7r[X^^^)} on ^ U {d} is the same 
as that of the Markov chain on ^ U {d} generated by Q^. By the assumption (A) we 
notice that < —X/q. Thus by using the strong Markov property we get 



m,0-\x,y)\ < -\\{Q^ + V^)-'\\^E^ 



(A.13) 



where is the first hitting time at y of the chain {tt{xP)}, and \\{Q^ + V^) ^||oo 
is the operator norm of {Q^ + ^)~^ acting on /°°(^ U {(?}), which satisfies the bound: 



\m+v^r'\\oo<q/\- 



(A.14) 



Since tP = for the Markov process starting at x, we see from ()A.13P and ()A.14p 
that 



m,o-\x,x)\<i/x. 



(A.15) 



Glauber and Kawasaki dynamics 



19 



This is the bound in the lemma for the diagonal components. In order to estimate 
|^(^, y)| for X ^ y, we define 



(A.16) 



If we let r the random variable which is exponentially distributed with parameter 1, 
then we have the identity: 



ulfHx) 



E 



-X/qr 



X + q 



%ix,y)+ %ix,z)ul^\z) 

%{x,y)+ %{x,z)u'y^\z) 
z&i-z7^y 



E 



1 .n,- 



{x,y). 



^ VA + g 

71=1 

Now by the definition of Q^{x, y) in (jA.lip we see that for ^ C 
n^(a;,y) = n^/(a;,y) whenever x,y € ^. 



(A.17) 



(A.18) 



Also, since d plays as a cemetery, once the process visits 5, it never comes out from 
it. Therefore, the nonzero contributions in the term (ng)"(x,y) come only from the 
random walk path of length n connecting x and y on the set ^ (avoiding the cemetery 
d). Obviously, such a number of paths increases as the set ^ increases. Therefore the 
last presentation of (IA.17P is bounded by T{x, y) where 



n=l 



\ + q 



-n 



(A.19) 



with n := n^;. Inserting this and (|A.14p into ()A.13p we get the bound in the lemma 
for off-diagonal components. Together with (jA.lSP we are done in the case that A is 
a real matrix. When ^ is a complex matrix, we write A = Ai + iA2, where Ai and 
A2 are real matrices. Let Ei and E2 be two copies of E. Then we have the bijection 



l\E) 9 / = /i + ^ /i © /2 G ll^,{E,) © ll,,{E2) = lL,{Ei U E2) 



(A.20) 



where ^reai(') ^e^ns the real Hilbert space. Under this map, A in l'^{E) is equivalent 
to the matrix 

V^2 Ai J 

acting on I'^^.^iiEi U E2). Similary, for any subset C E, the submatrix ^(^,^) acting 
on /^(^) is equivalent to the submatrix 



(A.22) 
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of A acting on Z^cail^i ^ ^2), where ^j, i = 1,2, are again the copies of ^. Notice that 
the enlarged real matrix A satisfies the conditions of the lemma with E being replaced 
by El U E2. Let A{C, C)-^ = C + iD, where C and D are real matrices. Then we can 
check that 

C = P^,A{tO^'P^, and D = -P^,AiC,0^'Pi„ (A.23) 

where P^^, i = 1,2, are the projections on Z^caiC^i ^ C2) onto /rcai(^i)) ^ = 1)2, respec- 
tively. So, < |C(x,y)| + \D{x,y)\, and by (|A.23p and applying the 
result for the real case we also arrive at the conclusion for the complex case. □ 
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